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BACKGROUND OF THE INVENTION 

I. Field of the Invention 

10 

The present invention relates to CDMA (Code Division 
Multiple Access) cellular telephone and wireless data 
communications with data rates up to multiple Tl (1.544 Mbps) and 
higher (>100 Mbps) , and to optical CDMA. Applications are 
15 mobile, point-to-point and satellite communication networks. More 
specifically the present invention relates to a new and novel 
approach to the design of waveforms and filters using 
mathematical formulations which generalize the Wavelet concept to 
communications and radar. 

20 



II. Description of the Related Art 

25 

Multi-resolution waveforms used for signaling and/or 
filters and which are addressed in this invention are defined to 
be waveforms of finite extent in time and frequency, with scale 

30 and shift properties of multi-resolution over the time-frequency 
(t-f) space. These waveforms can also be referred to as multi- 
scale waveforms, and include multi-rate filtering and the 
Wavelet as special cases. The enphasis will be on digital 
design and applications with the understanding that these multi- 

35 resolution waveforms are equally applicable to analog design and 
applications . 
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Background art consists of the collection of waveform and 
filtering design techniques which can be grouped into six broad 
categories. These categories are: CI) least squares (LS) design 
5 algorithms for filters and waveforms that design to 
specifications on their frequency response, C2) analytic 
filters and waveforms which are specified by a few free design 
parameters that can be sub-categorized into current applications 
and theorectical studies, C3) combinations of CI and C2 for 

10 greater flexibility in meeting communications and radar 
performance goals, C4) special design techniques to control the 
noise levels from intersymbol interference (ISI) and adjacent 
channel interference (ACI) in the presence of timing offsets for 
multiple channel applications, C5) Wavelet filter design using 

15 scaling functions (iterated filter banks) as the set of design 
coordinates or basis functions, C6) filter and waveform design 
techniques for non-linear channels and in particular for 
operation in the non-linear and saturation regions of a high 
power amplifier (HPA) such as a traveling wave tube (TWT) or a 

20 solid state amplifier, and C7) LS dynamic filters derived from 
discrete filtering and tracking algorithms that include adaptive 
equalization for communications, adaptive antenna filters, 
Wiener filters, Kalman filters, and stochastic optimization 
filters. 

25 

Category CI common examples of LS digital filter design are 
the eigenvalue algorithm in "A New Approach to Least-Squares FIR 
Filter Design and Applications Including Nyquist Filters" and the 
Remez-Exchange algorithm in "A Coirputer Program for Designing 

30 Optimum FIR Linear Phase Filters' - The eigenvalue algorithm is a 
direct LS minimization and the Remez-Exchange can be reformulated 
as an equivalent LS gradient problem through proper choice of the 
cost function. Both LS algorithms use the FIR (finite impulse 
response) digital sarrples as the set of design coordinates. LS 

35 design metrics are the error residuals in meeting their passband 
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and stopband ideal performance as shown in FIG- 3. Category C2 
common examples of analytical waveforms and filters for system 
applications are the analog Chebyshev, Elliptic, Butterworth, and 
the digital raised-cosine, and square-root raised-cosine. For 
5 theorectical studies, common examples are polyphase multirate 
filters, quadrature mirror filters (QMF) , and perfect 

reconstruction filters. Although these theorectical studies have 
yet to yield realizable useful filters for system applications, 
their importance for this invention lies in their identification 

10 and application of ideal performance metrics for filter 
designs. Category C3 common digital example is to start with 

the derivation of a Remez-Exchange FIR filter and then up-sample 
and filter with another bandwidth limiting filter. This results 
in an FIR over the desired frequency band that is larger than 

15 available with the Remez-Exchange algorithm and with sidelobes 
that now drop off with frequency compared to the flat sidelobes 
of the original Remez-Exchange FIR. A category C4 common example 
is to select the free parameters of the category C3 filter in 
order to minimize the signal to noise power ratio of the data 

20 symbol (SNR) losses from the ISI, ACI, and the non-ideal 
demodulation- A second common example is start with a truncated 
pulse whose length is short enough corrpared to the symbol 
repetition interval, to accommodate the timing offsets without 
significant impact on the ISI and ACI SNR losses. This shortened 

25 pulse can then be shaped in the frequency domain. 

Category C5 Wavelet filter design techniques discussed in 
the next section will serve as a useful reference in the 
disclosure of this invention. Category C6 common example is the 
30 Gaussian minimum shift keying (GMSK) waveform. This is a 
constant amplitude phase encoded Gaussian waveform which has no 
sidelobe re-growth through a non-linear or saturating HPA- 
Category C7 common exairples are the adaptive equalization 
filter for communication channels, the adaptive antenna filter, 
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and the Kalman filter for applications including target tracking 
and prediction as well as for equalization and adaptive antennas. 

Minimizing excess bandwidth in the waveform and filter 
5 design is a key goal in the application of the C1,...,C6 design 
techniques for communications and radar- Excess bandwidth is 
identified as the symbol a in the bandwidth-time product 
BTs=H-a where the two-sided available frequency band is B and the 
symbol repetition interval is Tg. Current performance capability 

10 is represented by the use of the square-root raised^cosine 

(sq-rt r-c) waveform with a=0.22 to 0-4 as shown in FIG. 8. The 
goal is to design a waveform with a=0 within the performance 
constraints of ISI, ACI, passband, sideband, and passband 
ripple. This goal of eliminating the excess bandwidth corresponds 

15 to the symbol rate equal to the available frequency band l/Tg = 
B/(H-a) = B for a=0. This symbol rate l/T^^B is well known to be 
the maximum possible rate for which orthogonality between 
symbols is maintained. A fundamental performance characteristic 
of the new waveform designs is the ability to eliminate excess 

20 bandwidth for many applications. 

Scope of this invention will include all of the waveform 
and filter categories with the exception of category C4 special 
design techniques and category C7 LS dynamic filters. Emphasis 

25 will be on the category C5 Wavelets to establish the background 
art since Wavelets are multi-resolution waveforms that can 
eliminate the excess bandwidth and have known design algorithms 
for FIR waveforms and filters. However, they do not have a design 
mechanism that allows direct control of the ISI, ACI, passband, 

30 sideband, and passband ripple. Category C2 theorectical studies 
also eliminate the excess bandwidth- However, they are not multi- 
resolution waveforms and do not have realizable FIR design 
algorithms. So the eirphasis in background art will be on 
Wavelets whose relevant properties we briefly review. 
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Wavelet background art relevant to this invention consists 
of the discrete Waveform equations and basic properties, 
application of Wavelets to cover a discrete digital time- 
5 frequency (t-f) signal space, and the design of Wavelets using 
the iterated filter construction. Wavelets are waveforms of 
finite extent in time (t) and frequency (f) over the t-f space, 
with multi-resolution, scaling, and translation properties. 
Wavelets over the analog and digital t-f spaces respectively are 
10 defined by equations (1) and (2) as per Daubechies' s "Ten 
Lectures on Wavelets", Philadelphia :_SIAM, 1992 



Continuous wavelet 

15 <^...(0=H""V(— ) (1) 

a 

Discrete wavelet 

i^..(«)=H"'V(— ) (2) 

a 

where the two index parameters "a,b" are the Wavelet dilation and 
20 translation respectively or equivalently are the scale and shift. 
The v|/ is the "mother" wavelet and is a real and symmetric 
localized function in the t-f space used to generate the doubly 

indexed Wavelet \|/ab- The scale factor "|a| '^^" has been chosen to 

keep the norm of the Wavelet invariant under the parameter change 
25 "a,b". Norm is the square root of the energy of the Wavelet 
response. The Wavelets \|/a,b and \|/ are localized functions in the 
t-f space which means that both their time and frequency lengths 
are bounded. The discrete Wavelet has the time "t" replaced by 
the equivalent digital sample number "n" assuming the waveform is 
30 uniformly sampled at "T" second intervals. 
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Wavelets in digital t-f space have an orthogonal basis that 
is obtained by restricting the choice of the parameters "a,b" to 
the values a=2~^ , b=qM2^ where ^p,q" are the, new scale and 
translation parameters and "M" is the spacing or repetition 
5 interval Ts=MT of the Wavelets (which from a communications 
viewpoint are symbols) at the same scale "p". Wavelets at "p,q" 
are related to the mother Wavelet by the equation 

10 t/p.,(n) -2-^ti^(2-Pn-qM) (3) 

where the mother Wavelet is a real and even function of the 
sample coordinates at dc (dc refers to the origin f=Q of the 
15 frequency f space) . The orthonormality property means that these 
Wavelets satisfy the orthogonality equation with a correlation 
value equal to ^^1". 

S<^P.q<^k.m =1 iff both p = kandq = m 

n 

= 0 otherwise (4) 

20 Wavelet representation of a digital t-f space starts with 

selecting an N sample time window of a uniform stream of 
digital samples at the rate of 1/T Hz (1/second) equivalent to a 
"T" second sampling interval. The N point or sample t-f space 
in FIG. 1 illustrates a Wavelet representation or "tiling" with 

25 Wavelets that are designed analytically or by an iterated filter 
construction. 

The t-f space in FIG. 1 is partitioned or covered or tiled 

by a set of Wavelet subspaces ^^,p = 0,l,...,W-l} 1 where N=2"'. Each 

30 Wavelet subspace Wp 1 at scale "p" consists of the set of Wavelet 
time translations {q = 0,1,..., N/2^*^-l} over this subspace. These 
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Wavelet subspaces are mutually orthogonal and the Wavelets within 
each subspace are mutually orthogonal with respect to the time 
translates. This N-point t-f space extends over the time 
interval from 0 to (N-l)T 2 and over the frequency interval from 
5 0 to (N-1) 3 in units of the normalized frequency fNT 4. 



The iterated filter bank in FIG. 2 is used to generate the 
Wavelets which cover the t-f space in FIG. 1. Each filter stage 
5 consists of a high pass filter (HPF) and a low pass filter 
10 (LPF) . Output 6 of the LPF is subsampled by 2 which is 
equivalent to decimation by 2 . This t-f space space is an N- 
dimensional complex vector metric space F 7. At stage m in the 
iterated filter bank, the remaining t-f space V^_^ 8 is 

partitioned into V^^j 9 and the Wavelet subspace Wm+i 10 . 

15 

Scaling functions and Wavelets at each stage of this filter 
bank satisfy the following equations 



20 ' (5) 

V^(«)=2-"2 5;g,(Z,(2n-q) Vq 



where (p is the scaling function, \|/ is the Wavelet, HPFp 

coefficients are {Z?^, Vq}, LPFp coefficients are {^g, Vq}, and 

the equations apply to the stages 0,l,...,m-l. Identifying the 
25 scale parameter and using the previous Wavelet formulations 
enable these equations to be rewritten for stages p=0, 1, ...,m-l 
as 

a 

(6) 
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For the application the HPFp and LPFp are quadrature mirror 
filters (QMF) with perfect reconstruction. This means they cover 
the siobspace Vp with flat responses over the subband frequency 
5 including the edges of the frequency subband, and the HPFp 
coefficients are the frequency translated coefficients for the 
LPFp: {gq=(-l)% , Vq}. 

Wavelet design using iterated filter bank starts with the 
10 selection of the scaling functions. Starting with a primative 
scaling function such as the one proposed by Daubechies-^y one can 
use the iterated filter construction given by equations (5) and 
(6) to derive successive approximations to a desired scaling 
function which has properties that have been designed into it by 
15 the selection of the filter coefficients {g^,Vq} at each level of 

iteration. The Wavelets can be derived from these scaling 
functions using the iterated filter construction or scaling 
equations (5) and (6) . 

20 Another use of the iterated filter construction is to 

design the scaling functions as Wavelets thereupon ending up with 
a larger set of Wavelets for multi-resolution analysis and 
synthesis as illustrated by Coifman' s Wavelets in "Wavelet 
analysis and signal processing" - 

25 



SUMMARY OF THE INVENTION 

30 

This invention for the design of multi-resolution Wavelet 
waveforms improves their performance for engineering and 
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scientific applications, and in particular for applications to 
communications and radar. Current practice is to 1) design 
traditional multi-resolution waveforms in the time domain using 
metrics which specify frequency performance and without 
5 consideration of Wavelet properties or to 2) design Wavelet 
multi-resolution waveforms using the time-domain iterated multi- 
resolution filtering approach with a set of scaling functions 
used to perform the filtering and without direct considerations 
of the frequency performance. These two separate design 

10 approaches yield fundamentally different waveforms. This 
invention provides a means to combine these two approaches to 
generate a new multi-resolution Wavelet waveform with the best 
properties of the traditional multi-resolution waveforms and the 
Wavelet multi-resolution waveform. 

15 

This invention introduces innovations for the design of the 
Wavelet waveform that 1) include the ISI (intersymbol 
interference), ACT (adjacent channel interference) , and QMF 
(quadrature mirror filter) requirements into the design algorithm 

20 along with non-linear modifications to the traditional passband 
and stopband frequency requirements which are used in the Remez- 
Exchange and eigenvalue LS algorithms, 2) select the frequency 
harmonics as the design coordinates, 3) use an appropriate 
subset of the available Fourier domain frequency harmonics as the 

25 design coordinates with the property that this subset of 
coordinates is a basis for the multi-resolution Wavelet waveform 
design, and 4) design the multi-resolution Wavelet waveform for 

no excess bandwidth a=0. 

30 These innovations support the development of the new 

waveforms 1) that are generalizations of Wavelets in the 
frequency domain using a means which makes them useful for 
communications and radar over their t-f space , 2 ) that can be 
designed to be an orthonormal (orthogonal and normalized) basis 
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or orthononnal set of coordinates over the available frequency 
bandwidth with the implicit property that the excess bandwidth 
vanishes a=0 to within the accuracies allowed by communication 
and radar design implementations, and 3) that have multi-scale 
5 properties which allow a single dc waveform design to be used to 
uniquely define the complete set of waveforms over the t-f space 
for multi-scale applications^ Specific design algorithm examples 
developed in this invention disclosure are the LS design 
algorithms that use eigenvalue and gradient search techniques in 
10 the frequency domain to find the best waveform design which 
minimizes the corresponding LS error residual cost function that 
is a weighted linear sum of the residual error metrics. 

Wavelet waveform performance calculated in this invcnoion 
15 invention disclosure from the application of these LS design 
algorithms, demonstrates the capability of this invention to 
provide a means to design waveforms that are improved over 
current practice. There are other algorithmic design algorithms 
which can be realized by modifications to this invention. ftn 
20 Two exanples given in this invention disclosure is — are the 
modification of these LS algorithms for application to constant 
amplitude bandwidth efficient (BEM) waveforms for communications 
and to synthetic aperture radar (SAR) waveforms . 

25 

BRIEF DESCRIPTION OF THE SEVERAL 
VIEWS OF THE DRAWINGS 

30 

The above-mentioned and other features, objects, design 
algorithms, and performance advantages of the present invention 

10 



will become more apparent from the detailed description set forth 
below when taken in conjunction with the drawings and performance 
data wherein like reference characters and numerals denote like 
elements, and in which: 

5 

FIG, 1 is a Wavelet N-point t-f space extending over the 
time interval (0, (N-l)T] and the frequency interval (0, (N- 
1)/NT]/ which is tiled or covered by a set of orthonormal 
Wavelets at the scales p=0, 1, ...,m-l . 

10 

FIG, 2 is a Wavelet iterated filter bank used to generate 
the set of Wavelets which tile or cover the t-f space in FIG. 1. 

FIG. 3 illustrates the power spectral density (PSD) of the 
15 new Wavelet waveform at dc, and the set of stopband and passband 
design requirements. 

FIG. 4 is the flow diagram of the LS metrics and cost 
functions and the final cost function used to find the optimal LS 
20 solution for the Wavelet waveform. 

FIG. 5 is the flow diagram for the eigenvalue recursive 
solution algorithm to find the optimal harmonic coordinates that 
define the LS solution for the Wavelet waveform. 

25 

FIG. 6 plots the dc PSD in dB units for the frequency 
response of the mother Wavelet at dc and the square^root raised- 
cosine (sq-rt r-c) waveforms with excess bandwidth parameter a= 
0.22, 0.40 vs. the normalized frequency in units of symbol rate 
30 I/T3. 

FIG. 7 plots the dc PSD in dB units vs. the normalized 
frequency fTg for the Wavelet waveform designed with 
modifications to the LS algorithms for application to constant 
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amplitude BEM communications, and for a Gaussian minimum shift 
keying (GMSK) waveform. 

FIG. 8 plots the amplitude of the dc radar ambiguity 
5 function vs. the normalized frequency fTp and normalized time tTc 
where the pulse time Tp and chip time Tc are both equal to the 
communications symbol interval Ts for this exartple, for the new 
waveform and an unweighted chirp waveform. 



DETAILED DESCRIPTION OF THE INVENTION 

15 New Wavelet waveforms in this invention disclosure are 

generalizations of Wavelets in t-f space which enable them to be 
useful for communications and radar applications- This 
generalization is accomplished by 1) the introduction of a 
frequency translation, by 2) changing the ortho normality 

20 condition in equation (4) to apply to waveforms within the same 
space {q} and over the scales {p} with the inclusion of the 
frequency translation, and 3) by their characterization and 
design in the Fourier domain. With frequency translation the 
analytical formulation of these new wavefomns as a function of 

25 the baseband or mother waveform centered at dc ( do refers to the 
origin f-0 of the froqucncy apace) b ecomes 

= 2""^ V^(2*'n - qM ) e^^^^-^^^^^^^"^ (7) 

30 

where fc(p,r) is the center frequency of the frequency translated 
dc waveform, at scale "p" and frequency index "r". The purpose 
of the frequency index ^^r" is to identify the center frequencies 
of the waveforms at the scale "p" in the t-f space. 

12 



These waveforms satisfy the complex orthonormality 

equations 



Ym 1^* =1 iff p = kandq = m and r = v 
5 =0 otherwise (8) 

where is complex conjugation, and are generalizations of the 

orthonormality equations for the analytical Wavelets in (4) . 

The new waveforms in equation (7) expand the Wavelet 
10 analytical formulation to include a frequency variable. 
Wavelets are functions of the scale and translation parameters 
''p,q"- The frequency variable together with the Fourier domain 
design are entirely new means for deriving these new waveforms 
as generalization of the traditional Wavelets 

15 

These waveforms are generalizations of Wavelets in the 
frequency domain using a means which makes them useful for 
communications and radar in the t-f space. The basis vectors for 
this metric space V consist of a subset of the admissible set 

20 of scaled and translated waveforms { V^p^g,r>^PjQj^ J derived from the 
dc waveform y/ as per equation (7) . An admissible waveform is 
any combination that covers f^=t-f space. We are interested in, 
the Fourier domain representation of the dc waveform y/ Ln V , 
and in particular in a subset of the discrete Fourier transform 

25 (DFT) harmonic coefficients over the Fourier domain which we 
intend to use as the design coordinates. Starting with the z- 
transform and continuous Fourier transform, the DFT harmonic 
coefficients are defined by the following equations 



30 
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DFT Harmonic Coefficients 



(9) 



= Hv^i^) z-transfonn 

n 

= Z<^(n) e"''*^ Fourier transform 

n 

= ZKn) W^r^ DFT harmonic coefficients for Vk 

n 

where 

k n 

^ = St^k sin((fi?/2- ;zk/N')Ar')/ Nsin(fi;/2 -;dc/lSP) 

k 

X [Harmonic interpolation for " k"] 

k 

= DFT frequency or harmonic coefficients such that = k/JfT 
where is the harmonic frequency corresponding to "k" 
N' = length of i//(n) 

wherein ^(cy) = is the power spectral density of the dc 

waveform. These equations define the frequency representation of 
10 the new waveforms in terms of the available set of harmonic 
coefficients, which set is considerable larger than required for 
most applications. 

The new waveforms are an orthonormal basis with no excess 
15 bandwidth, which properties are asymptoticly approached by the 
new waveforms to within design accuracies inherent in 
communications and radar. To demonstrate these properties we 
need to identify the structure of the dc waveform ±n V . We 
start with definitions for the parameters and coordinates in the 
20 following equation (10) . The waveforms derived from the dc 
waveform will be designed to be orthogonal over both time 
translates "MT" and frequency translates "1/MT" which 
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respectively correspond to the Wavelet symbol spacing Ts=MT and 
the adjacent channel spacing 1/Tg. This means the orthogonal 
spacing of the waveforms in V are at the time-frequency 
increments (MT, 1/MT) = (T3- l/T^) . In the interests of 
constructing the orthonormal multi-resolution waveforms to cover 

V it will be convenient to assume that M,L are powers of 2. We 

need the following definitions for the parameters and 
coordinates . 

Parameters and Coordinates (10) 

N' = Length of which is an even fixnction about the center and 
which spans an odd number of pomts or samples 
= ML + 1 where M, L are assumed to be even functions for 

convenience of this analysis 
= Number of pomts of i// 
M = Sampling interval for i// 

= Spacing of i// for orthogonality 
L = Length ofi// \n units of the sample interval M 

= Stretching of (// over L sample intervals 
n = n^ + njM 

= partitioning into an index n ^ over the sample 
length no =0,1,...,M -land an index nj =0,L..,L-1 
over the sample intervals 
k =ko + kjL 

- partitioning into an index ko over the harmonic 
frequencies ko =0,1,...,L-1 corresponding to the 
stretching and an mdex k^ = 0,1, ..,M over the harmonics 
frequencies corresponding to the admissible frequency 
slots for <^ 

The harmonic design coordinates are selected using the 
following observation. For most applications and in the following 
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development it is assumed that the waveforms are spectrally 
contained in the frequency interval 1/MT corresponding to the 
frequency spacing. This suggests the harmonic design coordinates 
be restricted to the subset of L harmonics {ko=0, 1, . . . ,L-1} 
covering this spacing. These L harmonics correspond to the 
stretching of the mother waveform over the L repetition 
intervals. 

Obviously, for some applications as will be demonstrated 
later, the spectral containment is spread out over several 
MT frequency increments whereupon one must increase the subset of 
design harmonics to possibly 2L, 3L or larger. 

The DFT equations for the dc waveform in (9) when rewritten 
in terms of the L harmonic design coordinates iWk^y } become: 

DFT equations for dc waveform 

^Jto ~ y^<^(n) W^**** harmonic design coordinates 

n 

y/{ri) ^(l/N^) ^V^k^ new waveform defined in terms (11) 

of the L harmonic design coordinates {^^/J^^,VkQ} 

It will now be shown that the use of these L harmonic design 
coordinates is sufficient under time translates to be a basis for 
the corresponding subspace of V which means these waveforms 
provide a complete set of coordinates to describe this subspace. 
We _thG will use the theorems of Karhunen-Loeve and Mercer and 
will limit the demonstration to the dc waveform for simplicity 
and without loss of generality. We start by considering the 
expansion of a random complex sequence {z(n), Vn} in a series of 
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waveform coordinates consisting of time translates of the mother 
waveform y/ . The sequence {z(n), Vn} is a zero-mean stationary 
random process which is orthonormal over the sample interval "M" 
and has a frequency spectrum which is flat and extends over the 

frequency range MT w hich is centered at baseband 

corresponding to a zero frequency- This means the {z(n), Vn} 
cover the subspace of V corresponding to the scale of the dc 
waveform if/ and its time translates {il/{n-qM) ^ {f/^{n\ . In 

addition, the V is now considered to be extended over a time 
interval which is relatively large compared to the N-dimensional 
t-f space in FIG. 1,2 to avoid end-effects on the analysis- We 
start by approximating the sequence {z(n), Vn} by the {z(w),Vn} 
where : 



where the complex coeflBcients {z^ }are 
derived from the original sequence using the 
orthogonality properties of i// 

n 

= 2 z(n) 

n 



The following equations prove that the coefficients {Zq, Vq} are 
orthonormal: 
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Z,Z^.* = Zz(An-qM)2*(An-q'M)ti^,(Aii)v^^. *(An) 



(13) 



= ^qq' Z|v^q(n)r siiice thc scquencc {z} is orthonormal 

n 

for qM time translates 
= {5_ , with normalization of the energy of y/ 



Equations (12) and (13) together prove the 



Karhunen- 



Loeve' s theorem which proves the following equation for the 
accuracy in approximating the stochastic sequence {2(n), Vn} by 
{z(w),Vn}. This accuracy is expressed by the expected "E(o)" 
squared error 'Mo)" in this approximation: 



We need to prove that the right hand side of this equation is 
zero which then proves that the approximating sequence is equal 
to the original sequence in the mean-square sense. In turn this 

proves that the new waveform coordinates { y/^^ Vq } are a basis for 
the original sequence {2(n), Vn} which is the goal. 

The right hand side of equation (14) when set equal to zero 
expresses Mercer's theorem so the goal is to prove Mercer's 
theorem. To do this we use the DFT of {// in equation (11) and the 
coordinates in (10) to evaluate the right hand side of equation 
(14) . We find 




(14) 
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where T = (l/LJ^W^'"^ 

= sin(;r Ako)/Zsin(7r Ak(j/L) 
= 1 for Ako =0 (15) 
= 0 otherwise 
This proves that 

1- Zk'^ -^'^i^f =^ ^'^o 



which proves that equation (14) reduces to 

5 

E\z{ri)-z{nf]^ 0 (16) 

Which^ as per the above_j_ proves that the set of multi-resolution 
waveforms is a basis. This proof easily generalizes to the 
10 multi-resolution waveforms at all of the scales {p} and time 
translates {q}, and to the expansion of the harmonic design 
coordinates over 2L, 3L, ... as required by the application. 

Next we demonstrate that these waveforms have multi-scale 

15 properties. that — we — will — dcmonotrato . ^We start with the 

observation that the Fourier domain design for the new waveforms 
provides a natural and easy way to derive the conplete set of 

waveforms {^p^,Vp,q} for the space V from the design of the dc 
waveform y/ by using the same invariant set of Fourier domain 
20 harmonic design coordinates {(//j^^^VIcq} derived for the dc multi- 
resolution waveform. This demonstration requires that we show 1) 
how the multi-scale transformations are implemented with the 
design in the Fourier domain^ and 2) how the waveform design 
remains invariant under scale changes. 

25 
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First consider the multi-scale transformation which derives 
the waveforms at the scale and shift parameters "p,q" from the dc 
waveform at scale ''p=0" and centered at the origin ^^q=0". We 
begin by extending the parameters and coordinates in equation 
(10) , to include both scaling and subsampling or decimation in a 
form that is equivalent to the iterated filter bank construction 
which is used to derive the current Wavelet waveform using the 
filter scaling functions. Starting with the coordinates at scale 
^'p=0" the parameters and coordinates at scale "p=p" are given by 
the equations: 



p-0 parameters and coordinates 

n = no + ni M (17) 

no = ao + ai 2 + ... + am-i 2^''^ 

= M points 
M = 2^" 

ni = bo + bi 2 + ... +bi-i 2^"^ 

= L points spaced at M sample intervals 

p=p parameters and coordinates 

2'P n(i2^) = scaled by "2"^" 

and subsampled or decimated by 2^:1 
wherein 2^:1 reads ""2^ to 1 decimation^^ 
= no(p) + ni(p)M 

2"P no(>l2^) = scaled by "2"^" 

and subsampled or decimated by 2^:1 
= no(p) 

= ap + ap-i2 + ... + ap+ni-i 2^^"*' 

= M points spaced at 2^ sample intervals 

2"P nx(>l<2^) = scaled by ^^2"^" 

and subsampled or decimated by 2^:1 
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= ni(p) 

ip + bp+i ^ T ... -r Up+i-i 
= L points spaced at M2^ sample intervals 



- bo + bo+i 2 + ... + bo+i-i 2^"^^"^ 



5 together with the observation that the sampling interval "T" is 
increased to "2'^T" under the scale change from "p-0" to "p=p" and 
subsampling or decimation from "1:1" to "2^:1". Combining these 
equations with the analytical formulation in (7) and the Fourier 
domain representation in (11) enables the waveforms at the 
10 parameters "p,q" to be written as a function of the Fourier 
domain harmonic design coordinates: 



= {2-"" /N') w;;.'^'*^-*^ e 



(18) 



*0 



for all admissible scale, translation, and frequency index 
parameters "p, q, r" . 

Next we need to demonstrate that the frequency domain 
20 design in (11) remains invariant for all parameter changes and in 
particular for all scale changes. This multi-scale property 
expresses the accordion behavior of the design in that the 
Wavelets at different scales are simply the stretched and 
compressed versions of the mother waveform with the appropriate 
25 frequency translation indices. This multi-scale invariancy means 
that the design for a M=16 channel filter bank remains the same 
for M=100 or M=10,000 channel filter banks, when the overlap L 
and the perf ojnnance goals remain constant. To demonstrate this 
invariant property across scales, we consider the multi- 
30 resolution waveform at scale "p" with the other parameters set 
equal to zero for convenience "q=0, r=0" and without loss of 
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generality. The Fourier domain frequency response can be 

evaluated starting with the original formulation in equation (7) : 



DFT at ^^p, q=0, r=0" (19) 

ko n(p) 



sin(;r (fN*2PT-ko)) 



_N'sin(;r (fNP2PT-ko)/N') 
= 2 V^ko [Hannonic mteq)olation for " "] 



*0 



This only differs from the harmonic representation in equation 
(9) in the restriction of the design coordinates to the subset of 
10 harmonic coefficients {t/^^jVkp} and the stretching of the time 

interval to "2PT" corresponding to the scale "p=p". The harmonic 
interpolation functions are observed to remain invariant over 
scale changes upon observing that the frequency scales as 
"f-'l/2PT" wherein ^^--^ is a proportionality operator^ w hich means 

15 the "frequency-time" product remains invariant with scale changes 
as per the fundamental property of the waveforms. This means the 
harmonic interpolation functions remain invariant with scale 
change and therefore the frequency response remains an invariant. 
This demonstrates the waveform design is an invariant across the 

20 waveform scales which means we only need a single design for all 
scales or resolutions of interest. 

LS design algorithms for new waveform will be described to 
illustrate the advantages the waveform has over current designs. 
25 The two LS algorithms described are the eigenvalue and the 
gradient search which respectively can be reduced to algorithms 
which are equivalent to current eigenvalue and Remez-exchange 
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waveform design algorithms for application to a uniform filter 
bank. We consider the t^f space which is spanned by a uniform 
polyphase filter bank consisting of M channels at the frequency 
spacing ^„=^_1/MT where T is the digital sampling interval, and 
5 the filter waveform FIR time response is stretched over L 
sampling time intervals Ts. This polyphase filter bank is ideally 
decimated which means the filter ouput sample rate I/T3 is equal 
to the channel-to-channel spacing T3=MT, equivalent to stating 
that there is no excess bandwidth a=0. The design for this 
10 topology is immediately applicable to an arbitrary set of multi- 
resolution filters through the scaling equation (18) which gives 
the design of the waveform at arbitrary scales in terms of the 
design of the dc waveform. 

15 For this polyphase filter bank used to construct the dc 

waveform or filter impulse response, the LS example design 
algorithms will use 5 metrics consisting of the 2 prior art 
passband and stopand metrics, and the 3 new metrics consisting 
of the ISI, ACI, and QMF, and solve the LS minimization problem 

20 using as design coordinates the subset of harmonic coordinates 
which are a basis. Since the two exaitple LS design algorithms 
only differ in the use of an eigenvalue LS optimization and the 
use of a gradient search LS optimization, the flow diagrams for 
the construction of the cost functions and the solution for the 

25 optimal waveform will be identical. However, there are 
differences in the construction of the cost functions from the 
respective metrics and in the iterative solution mathematics. 
Both LS solutions for the harmonic design coordinates minimize 
the weighted siom of the error residuals or cost functions from 

30 the 5 metrics. These design coordinates are the Fourier harmonics 

{^^j^^jVIcq} for {ko = 0,1,.„,L-1}. Resulting algorithms are easily 

extended to the applications requiring the design coordinates to 
cover 2L, 3L, ... hamonics. 
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Frequency domain design coordinates are related to the 
waveform time domain digital samples or coordinates as follows. 



Mappings of frequency to time (20) 

5 

Time domain design coordinates {\//(n),\fn} are real and 
symmetric and can be represented by the reduced set {ht(n), 
n=0,l,...,ML/2} 

ht(n) = i//(0) for n=0 
10 = 2i//{n) for n=l, 2, ...,ML/2 

= time domain design coordinates 

Frequency domain harmonic design coordinates { ^^^j^^ , Vk^ } are 

real and symmetric and can be represented by the reduced 
15 set {hf(k), k=0,l,,„.,L-l} 

hf(k) = l//^^ for k=ko=0 

= 2if/f^^ for k=ko=l, 2,...,L-1 

= frequency domain design coordinates 

20 Mapping of the frequency coordinates {hf(k), k=0, 1,..., L-1} 

into the time coordinates {ht(n), n=0, l,...,ML/2} is defined 
by the matrix transformation to within a scale factor 
ht = B hf 

where 

25 hf = (hf{0), ...,hf (ML/2) ) ^ transpose of column vector 

ht = (ht(0), ...,ht (ML/2) )^ transpose of column vector 
B = (ML/2 +1) X L matrix 

= [ Bkn ] matrix of row k and column n elements Bkn 
Bkn = 1 for n=l 

30 =2 cos(27ckn/ ML) otherwise 
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wherein the N'=ML+1 has been replaced by ML since a single end 
point has been added to the the FIR to make it symmetrical for 
ease of implementation for the example Wavelet being considered_j_ 
with a sample at the mid-point that makes the number of sartples 
5 N' an odd number. 

Passband and stopband metrics and cost functions are 
derived with the aid of FIG. 3 which defines the pe — PSD 
parameters of interest for the passband and stopband of the PSD 

10 for communications applications. Requirements for radar 

applications include these listed for communications. Referring 
to FIG. 3 the passband 11 of the waveform PSD is centered at dc 
(f=0) since we are designing the dc or baseband waveform, and 
extends over the frequency range ©p extending from -(Dp/2 to +a)p/2 

15 12 in units of the radian frequency variable (0=27cfT 13. The 
frequency space extends over the range of f=-l/2T to f=+l/2T 
which is the frequency range in FIG 1 translated by -1/2T so 
that the dc waveform is at the center of the frequency band. 
Quality of the PSD over the passband is expressed by the passband 

20 ripple 14 . Stopband 15 starts at the edge 16 of the 
passbands of the adjacent channels +/-(0a/2 16 and extends to 
the edge of the frequency band (o=+/-7t 17 respectively. 
Stopband attenuation 18 at +/-(Da/2 measures the PSD isolation 
between the edge of the passband for the dc waveform and the 

25 start of the passband for the adjacent channels centered at +/-ca3 
19. Rolloff 20 of the stopband is required to mitigate the 
spillover of the channels other than the adjacent channels, onto 
the dc channel. Deadband or transition band 21 is the interval 
between the passbands of contiguous channels, and is illustrated 

30 in FIG. 3 by the interval from G)p/2 to a)a/2 between the dc 
channel and adjacent channel at ©a- Waveform sample rate ©3 22 
is the waveform repetition rate. For the LS example algorithms, 
the waveform sample rate is equal to the channel-to-channel 
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spacing for zero excess bandwidth. Therefore, 1/Ts== (03/27tT=l/MT 
which can be solved to give 0)s=27c/M for the radian frequency 
sampling rate of the filter bank which is identical to the 
waveform repetition rate. 

We start by rewriting the DFT equations for the dc wavefom 
in (11) as a function of the {hf(k), k=0, 1, L-1 } 

~ ^t/(n)cos(nty) 

n 

= B hf using (21) and the definition of the vector "c" (21) 
c = ( 1, cos{<d\ cos( (ML/2) J transpose of column vector 

which is equivalent to the equation for y/{Q>) in (9) expressed as 
a linear function of the {hf(k), k=0, 1, L-1 } . An ideal "c" 
vector "Cj:" will be introduced for the passband and the stopband 
in FIG. 3 in order to identify the error residual Slf/(0) at the 
frequency "^y" in meeting the ideal passband and stopband 
requirements. The ideal PSD is flat and equal to "1" for the 
passband, and equal to "0" for the stopband. We find 

Error residuals for passband and stopband (22) 

Cr = (1,1,.», 1)*^ passband ideal "c" 
= (0,0,..., 0)*" stopband ideal ^^c" 
Sc = Cr - c error vector 

Si//i0) = &'Bhf 

= Residual error relative to the ideal spectrum at " " 
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The LS metric for the passband and stopband can now be 
constructed as follows for the eigenvalue and the LS optimization 
(or equivalently, the LS algorithm) design algorithms 



Passband and stopband metrics (23) 



J (band) = — ^ — f|(J(^(fi?)| do Eigenvalue 
^ Wj^'R hj^ Eigeinvalue 

- ¥y IP 



LS 



where 

band = [0, (Op ) passband 
~ (©3 , 7c] stopband 

R = \(&Sc5c'B)d0 

= L X L matrix 
5\|/ = { 5\|/(o)i) 5\|/((Du) )^ 

= vector of error residuals at the 

frequencies ©i , Ou across the band 

= norm or length of the vector (o) and which 

includes a cost function for the errors of 
the individual components 



where it is observed that the eigenvalue approach requires that 
the LS metrics be given as quadratic forms in the design 
coordinates {hf(k), k=0, 1,.„,L-1} whereas with the LS approach it 
is sufficient to give use norm-squared the LS metrics^^ qo' vector 
norm. 
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QMF metrics express the requirements on the deadband that 
the PSD' s from the contiguous channels in FIG. 3 add to unity 
across the deadband [cOp , ©s] in order that the filters be QMF 
5 filters. By suitable modification of the error vector 8c, the 
previous construction of the passband and stopband metrics can be 
modified to apply to the deadband. We find 
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Deadband metrics 



(24) 



15 



J (deadband) = ^ fl^^V^C^y)!^ do) Eigenvalue 



deadband 



hi' 



where 



5c 



= Cr - C((0) - C{7C/M - <0) 



(24) continued 
Eigenvalue 

LS 



20 

where c(o)=c as defined in (22) and (23), and c(7c/M - ay) = c at 
the offset frequency "tc/M - o" corresponding to the overlap of 
the contiguous filters over the deadband. 

25 Orthonormality metrics measure how close we are able to 

designing the set of waveforms to be orthonormal over the t-f 
space, with the closeness given by the ISI and the ACT. ISI and 
ACI errors are fundamentally caused by different mechanisms and 
therefore have separate metrics and weights to specify their 

30 relative importance to the overall sum of the LS metrics. ISI 
is a measure of the non-orthogonality between the stream of 
wavefoms within a channel as per the construction in FIG. 3. 
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On the other hand, ACI is a measure of the non-orthogonality 
between the waveform within a channel and the other waveforms in 
adjacent channels. This means the stopband performance metric 
has a significant impact on the ACI due to the sharp rolloff in 
frequency of the adjacent channel, and the ACI metric is then a 
measure of the residual non-orthogonality due to the inability of 
the stopband rolloff in frequency from completely eliminating the 
ACI errors. 

We assume that the received waveform is identical to the 
filter waveforms and is transmitted at the filter output sample 
intervals equal to MT seconds. The second assumption means we 
are assuming the receiver is synchronized with the received 
signal. Since there is no information lost by sampling 
asynchronously with the received waveform, we are free to make 
this synchronization assumption without loss of generality. ISI 
metrics are derived in the following set of equations. 

ISI metrics (25) 
Mapping of hf into \|/ 

\|/ = (\|/(-ML/2), ... , \|/(ML/2))^ transpose of column vector 
= H hf 

H = (ML+1) X L matrix of elements Hjcn 

Hkn = 1 for n=l 

= 0.5 cos (27Ckn/ML-l-l) otherwise 



Offset matrix A 
A = L X (ML+1) 
Akn - 0 

= [0 0 ... \|/(-ML/2) 



matrix of elements Akn 
for row vector k=l 

... \|/( (L-k)M+l) ] row vector k 
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12 kM LM+1 column elements 

ISI error vector 5e 

6E = L X 1 coliamn vector 



= A H h 



ISI metric 

J (ISI) = 8E^ 6e Eigenvalue 

= Non-linear quadratic function of hf 

10 = \\SE\\' LS 

ACI metrics are derived using the ISI metric equations with the 
following modifications • 

15 ACI metrics (26) 

Mapping of hf into \|f is the same as developed for ISI 
Offset matrix A elements are changed as follows to apply to 
Channel 1: 

20 Akn - 0 for row vector k=l 

- [0 0 ... \|/(-ML/2)PF^ ... \|/((L-k)M+l) W^^~'^] 
1 2 kM+1 LM+1 

which means the ACI error vector 5E is 
5e = L X 1 column vector 
25 = A H hf 

ACI metric for the two contiguous channels 
J (ISI) = 2 5E^ 5E Eigenvalue 

- Non-linear quadratic function of hf 

30 =2 pE\\' LS 
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where the factor "2" takes into account there are two contiguous 
channels or one on either side of the reference channel 0 in FIG. 
3. Because of the fast rolloff of the frequency spectrum the 
addition of more channels into the ACI metric is not considered 
5 necessary, although the functional form of the ACI metric in (26) 
allows an obvious extension to any number of adjacent channels 
which could contribute to the ACI. 



10 of the LS metrics derived in (23) , (24) , (25) , (26) . The LS 

algorithms minimize J by selecting the optimal set of frequency 
coordinates {hf(k),Vk} for the selected set of parameters used 
to specify the characteristics of the dc waveform, frequency 
design coordinates, LS metrics, and weights. Cost function and 
15 optimization techniques are given by the equations 



Cost function J for the LS algorithms is the weighted sum 



Cost function J 



(27) 



J 



^w{metrics) J(metrics) 



X 
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weighted sum of the LS metrics J (metrics) 



where 



metrics = passband, stopband, deadband, ISI, ACI 
{w (metrics ), Vmeir/C5} = set of weights 




= 1 normalization 



X 
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Optimization goal 

Goal: minimize J with respect to the selection of 



the {hf(k),Vk} 
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Optimization algorithms 

Two algorithms are the Eigenvalue and the LS optimization 



where the eigenvalue optimization algorithm uses the non- 
linear quadratic formulations of the LS metrics and the LS 
optimization algorithm uses the norm formulations for the 
LS metrics. 

5 

FIG. 4 is a summary of the LS metrics and the construction 
of the cost function J LS algorithms which use quadratic error 
metrics and LS algorithms which use norm-squared metrics. Design 

10 parameters 23 are the input and output design parameters. 
Input parameters are the number of polyphase channels M or 
equivalently the number of digital samples at spacing T over the 
symbol interval Ts=MT/ the length of the FIR time response for 
the waveform in units of L which are the number of digital 

15 samples per waveform repetition interval Tg so that the total 
number of digital samples for the symmetric FIR time response is 
equal to N'=ML+1, number of DFT samples per FIR length n_fft for 
implementation of the LS algorithms, passband radian frequency 
(Op, stopband radiam frequency ©ar waveform repetition rate in 

20 radian frequency Og, selection of the set of design coordinates 
{hf} to be used in the optimization, and the metric weights 
{w (metrics) } . Output parameters are the set of harmonic design 
coordinates {hf} that minimize J. Band metrics 24 are the 
passband, stopband, and deadband metrics defined in equations 

25 (23) , (23) , (24) respectively. Interference metrics 25 are the 
ISI and ACI metrics defined in equations (25) and (26) 
respectively. LS cost function J 26 is the weighted linear sum 
of the metrics defined for the band 24 and the interference 25 
as defined in equation (27) . 

30 

FIG. 5 is a flow diagram of the LS recursive solution 
algorithm using quadratic error metrics and norm-square metrics. 
There are two loops with the topology constructed so that the 
outer loop 27 is an iteration over the set of metric weights 
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{w (metrics) } , and the inner nested loop 28 is the recursive or 
iterative^ LS solution to find the optimal {hf} for the given 
design parameters and weights where each step refers to an 
iteration step in the solution. A recursive LS solution is 
5 required to be able to solve the highly non-linear interference 
metrics as well as the band metrics when one chooses to use non- 
linear techniques to construct these metrics. 

Consider the inner nested loop 28. The recursive LS 

10 solution starts with the initial step i=0 29 which begins with 
the selection of the input design parameters and weights and the 
selection of an initial set of values for the {hf} 30. Next the 
band metrics are calculated 31. Since this is the initial step 
i=0 32 the cost function J is restricted to the linear band 

15 metrics in order to fine the linear approximation to {hf} 33 to 
initialize the non-linear solution starting with step i=l 34 
wherein the highly non-linear interference metrics are included 
in J. Following the signal flow, for step i=l the band metrics 
31 and interference metrics 35 are calculated, weighted, and 

20 summed to form the cost function J 36. An iterative solution 
algorithm 37 finds the best approximation in step i=l to the 
{hf} which minimizes J using the approximation to {hf} from the 
previous step i=0 to linearize the search algorithm coefficients 
for the eigenvalue and LS step i=l iteration. If there is no 

25 convergence to the correct solution for {hf} this recursive 
sequence of calculations is repeated for step i=2 34. This 
recursive solution technique is repeated for subsequent steps 
until there is convergence 38 whereupon one exits this inner 
loop. 

30 

Consider the outer loop. After exiting the inner loop 38 
the solution for {hf} is tested to see if it meets the 
performance goals 39. If not, a new set of metric weights is 
selected and the inner loop is initialized i=0 40 and the 
35 inner loop is used to find the next solution for {hf}. This 



process continues until the performance goals are met or 
Qdoquatcd adequately approximated 39 whereupon the algorithm 
is exited with the final solution set {hf} 41. 



5 Applications of this new invention to both communications 

and radar will be given using these example algorithms and other 
algorithms supported by this invention . These new waveforms are 
considered for the applications: l)to replace the square-root 
raised^cosine waveform (sq-root rc) which is extensively used for 
10 the third generation (3G) CDMA communications, 2) to replace the 
Gaussian minimum shift keying (OMSK) waveform for constant 
amplitude bandwidth efficient (BEM) applications, and 3) as a 
candidate waveform for synthetic aperture radar (SAR) and real 
aperture radar (RAR) applications. 

15 

CDMA communications application for the current and the new 
3G CDMA considers a waveform designed with this new invention as 
a possible replacement for the sq-rt r-c waveform with bandwidth 
expansion parameter a = 0.22 to a = 0.40. This notation means that 

20 for a '022 the spectral efficiency is (symbol rate /bandwidth) = 
l/l+Of =1/1.22 = 0.82 = 82%. A basic advantage of the waveform 
is the potential for a symbol rate increase within the same 
bandwidth with an increase in the spectral efficiency to « 100% 
depending on the application and operational constraints. The dc 

25 power spectral density or power spectrum (PSD) of the waveform is 
compared to the PSD for the sq-rt r-c in FIG. 86. Plotted are the 
measured PSD in dB 42 versus the frequency offset from dc 
expressed in units of the symbol rate 43 Plotted against the 
normalized frequency offset are the dc PSD for the new waveform 

30 44, the sq-rt r-c with a=0.22 45, and the sq-rt r-c with 
a=0.40 46. It is observed that the PSD for the new waveform 
rolls off faster that that for the sq-rt r-c which means that the 
new mr waveform will support an increased symbol rate for a given 
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available frequency band while satisfying the inherent 
requirements for low ISI and MAI (multiple access interference) . 

Constant amplitude BEM application of the new waveform 

5 indicates that it is a viable candidate for replacing the current 
preferred modulation waveform which is the GMS K designed to mimic 
a Gaussian waveform by suitable modification of the constraints 
on the error metrics . The GMSK finds applications for 
transmitters which operate their HPA(s) aitplifiers in a 

10 saturation mode in order to maximize their radiated power from 
the HPA(s), and which require a BEM PSD to avoid excessive 
spreading of the transmitted power. Simulation data for the dc 
PSD is plotted in FIG. -9-7 for the new waveform BEM and the GMSK. 
Plotted are the measured dc PSD in dB 47 versus the frequency 

15 offset from dc expressed in units of the bit rate 48. Plotted 
against the normalized frequency offset are the dc PSD for the 
new waveform BEM 49 and the GMSK 50, for a length parameter 
L=10 where L is the length of the phase pulse in terms of the 
phase pulse repetition rate. The significance of this example 

20 data is that the new waveform has the potential to be designed 
to offer a PSD which is less spread out than the current GMSK, 
and therefore an improved BEM waveform. 

Real aperture radar (RAR) and synthetic aperture radar 
25 (SAR) applications of the wavefor m designed to mimic a Gaussian 
waveform by suitable modification of the constraints on the error 
metrics indicate that it is a viable candidate to replace the 
current chirp waveforms for wideband signal transmission, when 
combined with pseudo-random phase codes. Results of the 
30 simulation for the new waveform and an unweighted frequency 
chirp waveform are given in FIG. -iOS. Plotted are the ambiguity 
function for the new waveform 51 and the unweighted frequency 
chirp waveform 52. The dc 2-dimensional radar ambiguity 
function 53 is plotted as a function of the frequency offset in 
35 units of fTp and the time offset in units of t/Tc where Tp is 

35 



the phase-coded radar pulse length or length of the phase code 
and Tc is the phase code chip length. The chip length is 
identical to the waveform repetition interval Tg so that Tc = Tg. 
It is observed that the new waveform has the potential for 
significant improvements in the ambiguity function and by 
implication in the performance. 

Preferred embodiments in the previous description are 
provided to enable any person skilled in the art to make or use 
the present invention. The various modifications to these 
embodiments will be readily apparent to those skilled in the art, 
and the generic principles defined herein may be applied to other 
embodiments without the use of the inventive faculty. Thus, the 
present invention is not intended to be limited to the 
embodiments shown herein andis to be accorded the wider scope 
consistent with the principles and novel features disclosed 
herein. 
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5 ABSTRACT OF THE DISCLOSURE 

A method for designing Wavelets for communications and 
radar which combines requirements for Wavelets and finite impulse 
response FIR filters including no excess bandwidth, linear 

10 performance metrics for passband, stopband, quadrature mirror 
filter QMF properties, intersymbol interference, and adjacent 
channel interf erence, polystatic filter design requirements, and 
non-linear metrics for bandwidth efficient modulation BEM and 
synthetic aperture radar SAR. Demonstrated linear design 

15 methodology finds the best design coordinates to minimize the 
weighted sum of the contributing least-squares LS error metrics 
for the respective performance requirements. Design coordinates 
are mapped into the optimum FIR symbol time response. Harmonic 
design coordinates provide multi-resolution properties and enable 

20 a single design to generate Wavelets for arbitrary parameters 
which include dilation, down- samp ling, up-sampling, time 
translation, frequency translation, sample rate, symbol rate, 
symbol length, and set of design harmonics. Non-linear 
applications introduce additional constraints. Performance 

25 examples are linear communications, BEM, and SAR. 
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